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Abstract 

We consider the heat transfer problem associated with a periodic array of extended 
surfaces (fins) subjected to convection heat transfer with a uniform heat transfer coeffi¬ 
cient. Our analysis differs from the classical approach as (i) we consider two-dimensional 
heat conduction and (ii) the base of the fin is included in the heat transfer process. The 
problem is modeled as an arbitrary two-dimensional channel whose upper surface is flat 
and isothermal, while the lower surface has a periodic array of extensions/fins which are 
subjected to heat convection with a uniform heat transfer coefficient. Using the gener¬ 
alized Schwarz-Christoffel transformation the domain is mapped onto a straight channel 
where the heat conduction problem is solved using the boundary element method. The 
boundary element solution is subsequently used to pose a shape optimization problem, 
i.e. an inverse problem, where the objective function is the normalized Shape Factor and 
the variables of the optimization are the parameters of the Schwarz-Christoffel transfor¬ 
mation. Numerical optimization suggests that the optimum fin is infinitely thin and that 
there exists a critical Biot number that characterizes whether the addition of the fin would 
result in an enhancement of heat transfer. The existence of a critical Biot number was 
investigated for the case of rectangular fins. It is concluded that a rectangular fin is 
effective if its thickness is less than 1.64 k/h, where the h is the heat transfer coeffi¬ 
cient and k is the thermal conductivity. This result is independent of both the thickness 
of the base and the length of the fin. 
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1 Introduction 

Assuming isothermal boundary conditions, any extension from a surface would result in a 
reduction in heat transfer rate. This statement can be deduced from geometrical inclusion 
theorems [T] . It implies that the heat transfer rate across a two-dimensional channel, bounded 
by isothermal surfaces, is higher than that across a similar channel whose surfaces are ex¬ 
tended, i.e. by adding fins [2], assuming that the temperature difference between the surfaces 
remains the same. Hence, within the approximation of isothermal conditions and constant 
temperature difference, the addition of fins would simply reduce the heat transfer rate! Be¬ 
cause isothermal conditions can be realized in the limit of infinite Biot number, the following 
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Nomenclature 

Bi 

Biot number Bi = hL/k 

H 

dimensionless height of the slab/channel (dimensionless) 

Hb 

thickness of the base (dimensionless) 

Hf = H-Hb 

length of fin/extended surface (dimensionless) 

h 

height of the slab/channel in the transformed domain (dimensionless) 

h 

heat transfer coefficient (W/(m^ K)) 

i 

imaginary unit 

k 

thermal conductivity (W/(m K)) 

p 

arc-length of the periodic geometry (dimensionless) 

L 

period of the geometry (m) 

T 

temperature (dimensionless) 

x,y 

coordinates of the physical plane (dimensionless) 

z = X + iy 

complex coordinate of the physical plane (dimensionless) 

Zi 

vertices in the physical plane 

Greek symbols 

at 

turning angles 

w 

complex coordinate of the transformed domain (dimensionless) 

Wi 

image of Zi vertices in the transformed domain 

WN 

period of the geometry in the transformed domain 

Subscripts 

i 

related to the Tth vertex 

Diacritic 

A 

the variable is normalized with wn 


question is raised: When is an extended surface (fin) effective? This question will be 
elucidated in this paper, where we show that the limit of infinite Biot number is singular 
[H [4] , and a critical thickness exists [5] . 

This paper considers the fundamental problem of finding the optimum shape of an ex¬ 
tended surface/fin such that the heat transfer rate is maximized. In particular, we consider 
the inverse design problem associated with two-dimensional (2D) heat conduction in a finite 
2D periodic channel/slab with a flat isothermal upper boundary and a periodic lower bound¬ 
ary which is subjected to convection with a constant heat transfer coefficient, i.e. convection 
heat transfer is only considered to the extent that it provides a boundary condition for the 
conduction problem. The heat flux is proportional to the temperature difference between the 
surface and the far field [2] i.e. 


k dT/dn = h (Too - Tsurface)- 

Isothermal conditions can be realized in the limit of strong convection, i.e. large Biot number. 
Here, we should point out that isothermal conditions have been widely used in the heat transfer 
analysis of extended surfaces due to the very little information available on the coupling 
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between fin conduction and fluid convection, and the weak dependence of the heat transfer 
coefficient on the temperature difference between the base of the fin and its tip ([7]; [6], 4.36; 
i, 4.5). 

As we have mentioned, the optimization problem is an inverse design problem in the 
sense that the objective is to find the geometry that maximizes the heat transfer rate as 
oppose to the classical/direct problem, where the objective is to find the heat transfer rate 
associated with a given geometry. The objective function is the Shape Factor mi, i.e. the 
total heat transfer rate, and the variable of the optimization is the shape of the pipe which is 
parameterized though the parameters of the generalized Schwarz-Christoffel transformation. 
Hence, using Geometry Parametrization [El Emails], the Shape Optimization problem is 
posed as a nonlinear programming problem (constrained nonlinear optimization [Hj), which 
is solved numerically m- 

For regular, symmetric, isothermal, doubly-periodic walls, the heat conduction problem in 
a semi-infinite domain was addressed by Fyrillas &: Pozrikidis m using both boundary inte¬ 
gral and asymptotic methods. For 2D periodic channels/slabs, the problem has been addressed 
by Brady & Pozrikidis [E] where the authors considered the heat conduction problem asso¬ 
ciated with irregular isothermal periodic surfaces, using the generalized Schwarz-Christoffel 
transformation developed by Davis [E], Floryan [20j and Floryan & Zemach |21j . It was 
concluded that for regular geometries the shape plays a more important role in the total 
transport rate rather than the total arc-length while, for self-similar irregularities, the height 
of the roughness is the significant factor. These conclusions lead naturally to the follow¬ 
ing question: Given the arc-length and the period of a periodic surface/curve, what is the 
geometry that maximizes the overall transport rate; This Shape Optimization problem was 
addressed by Leontiou, Kotsonis &: Fyrillas [T8] . 

From an engineering perspective, knowledge of a surface geometry that maximizes the 
transport rate offers opportunities for new designs that exhibit enhanced characteristics and 
properties. For example, the problem of transport across an uneven surface, described in 
the preceding paragraph, is relevant to a variety of engineering applications involving heat 
transfer across rough and irregular boundaries, such as the surface of a circuit board in 
microelectronics [221 123] I25j . In general, heat transfer in slab-like configurations is of 

interest to problems associated with Heat Transport from Extended Surfaces (Fins) |8] and 
inverted high conductivity fins/inserts [26]. Published work in these thematic areas |27 ( I28 ( 
EniiMj suggests that there is potential for significant improvements if one considers a two- 
dimensional (2-D) heat conduction model as suggested by Aziz [3T] . 

In Section §2 we address the heat transfer problem associated with a periodic array of 
periodic extensions/fins of uniform convection heat transfer coefficient. In Section §3 we 
address the shape optimization problems associated with the optimum shape of the fins such 
that the heat transfer rate is maximized. In particular, in Section §3.2 we investigate under 
what conditions a rectangular fin enhances the heat transfer rate. 

2 Shape factor of a periodic array of extended surfaces (fins) 

In this section we consider 2-D heat conduction in a finite slab. The geometry of the slab is 
periodic in the horizontal direction and bounded in the vertical direction by an isothermal (Tq) 
flat surface at the top, while the bottom surface is subjected to convection with a uniform 
convection heat transfer coefficient (/)) and a constant far-held temperature T^o [2]- The 
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bottom surface is not flat, rather a periodic array of extensions is present (extended surfaces, 
fins), in order to enhance the heat transfer rate [21 [8]. Continuity of the heat flux at the 
bottom surface implies that the heat conduction rate, k = k dT/dn, must be equal to 

the heat convection rate h (T — Too), where k is the thermal conductivity. 

We non-dimensionalize lengths with the distance between the fins (period T), and the 
temperature by subtracting To and dividing by the temperature difference Too — Tq. The 
dimensional analysis leads to the following definition for the Biot number, Bi = L h/k. The 
domain and the dimensionless parameters associated with the problem are clearly indicated 
in Fig. [H 



Figure 1: Schematic representation of the problem in the physical domain. All variables are 
non-dimensional; lengths are non-dimensionalized with the distance between the periodic fins, 
i.e. period L. The dimensionless thickness of the base of the fin is Hi,, and the length of the 
is Hf = H — Hb. The non-dimensional temperatures are T = 0 at the top boundary and 
Too = 1 at the far field. 


2.1 Conformal transformation of the physical domain into a straight chan¬ 
nel 


To find the heat transfer rate of such a periodic slab we first transform it into a straight 
channel [3211331 (Ml 1181 ITlISSj. The relevant transformation, the generalized Schwarz-Christoffel 
transformation applicable for periodic channels, has been developed by Davis and Floryan 
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where the inner product identifies the number of elements (N) of the discretized lower bound¬ 
ary (fin), and the infinite outer product the periodic nature of the domain. We normalize 
lengths in the complex domain with Wj^, i.e. the period in the transformed domain. Hence, 
in the above transformation, wjs are the normalized images of the zjs vertices, a represents 
the A^-|-1-tuple 00 , 01 , 0 : 2 ,... ,aj^ which are equal to the turning angles multiplied by tt (the 
angles are taken to be positive for a clockwise rotation, and oq and Oj^ are defined with 
respect to the x-axis), T is a complex constant, and h is the normalized height of the chan¬ 
nel in the transformed domain (without loss of generality we assume that h = H). For the 
configurations we will consider i? is a real number and can be obtained by requiring that the 
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upper wall of the physical plane, i.e. the line z = iH, transforms to w = ih: 


Im 


\ih: 


a 


= H. 


(2) 


In addition, in view of the geometry, we must have 

N 

j=0 

Given the domain, the parameters of the transformation m can be calculated by solving 
a system of non-linear equations [iai2o]. However, in this work, we pose an optimization 
problem where the lower boundary is the variable of the optimization. In particular, we look 
for the optimum shape of the extended surface such that the heat transport rate is maximized. 
We pose the problem with respect to the parameters ct’s of the Schwarz-Christoffel transfor¬ 
mation. Essentially, we parameterize the lower boundary [ini El da 13] with respect to the 
parameters q’s, which are the variables of the optimization procedure, while the objective 
function is the heat transfer rate (the Shape Factor). An expression for the Shape Factor is 
obtained in the following section using the boundary element method. 


2.2 Shape Factor of an extended surface (fin) 

In view of the conformal transformation and the boundary condition on the lower surface we 
can obtain the following expressions for the Shape Factor (S) [2] associated with an extended 
surface of unit span: 


dT ^ 

S = -l ^[C,V = 0]<1C = B^ 


i-r[e,r) = o: 


drc 


d6 


(4) 


W=^ 


We also define the fin effectiveness (ej) as the ratio of the Shape Factor, as defined above, to 
the Shape Factor of the base without fins, i.e. Sb = Bijil -|- Bi Hb): 


£/ = 




(5) 


Hence, the definition of the Shape Factor and the fin effectiveness associated with a fin departs 
from the classical definitions [2], as it includes the area of the base not covered by the fin. 
Based on the above definition, an addition of an extended surface or fin would improve the 
heat transfer rate of the base, if its effectiveness is greater than one. 

The temperature along the lower surface, which includes the fin, can be obtained by 
applying the boundary element method [Ml EH El ESI ESI El E]. An appropriate Green’s 
function is that associated with a periodic array of sources of period 1 located along an 
insulated lower surface and a Dirichlet boundary condition along the top surface as described 

in [Ml EH 11 ES]: 


1 °° 1 
= 0] = — ( h -\— — tanh 

^ \ 'TT * ^ m 
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Figure 2: Mapping of the physical domain onto a straight channel using the generalized 
Schwarz-Christoffel transformation (equation [T]). Lengths in the physical and complex domain 
are non-dimensionalized with L and respectively. 


The mathematical domain along with the boundary conditions is shown in Fig. [3l Apply¬ 
ing the boundary element method and the boundary conditions we obtain a Fredholm Integral 
equations of the second kind for the temperature along the lower surface, 


T[i] = Bi 



G 


i' - i 


dw 


(r[|] -1) d|. 

w=^ 


(7) 


The integral equation is solved numerically using the collocation boundary element method 
IMlETliaEHlIMlIMlElISllIl!- It is important to note that in the above formulation the 
Generalized Schwarz-Christoffel transformation appears in the boundary element formulation 
similar to [35] . Another complication is that unlike [35| , the singularities of the transformation 
also appear in the boundary element formulation, which are addressed using Gauss-Jacobi 
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T = 0 



Figure 3: Schematic representation of the model problem along with boundary conditions. 
The lower boundary, which represents the extended surface, is subjected to convection heat 
transfer. Note that the conformal transformation has introduced an artificial heat transfer 
coefficient that reflects the transformation of the geometry. 


quadrature as explained in the following section, where we pose the Shape Optimization 
problem. The objective function is the fin effectiveness (e/), i.e. the normalized Shape 
Factor, and the variable of the optimization is the shape of the lower surface (fin); the latter 
is characterized through the parameters of the generalized Schwarz-Christoffel transformation. 

3 Optimum extended surfaces 

The formulation of the Shape Optimization problem follows along the same lines as the prob¬ 
lems formulated by Fyrillas and Leontiou & Fyrillas [32l |33l [Ml ISl IH ES]. The objective 
function is the fin effectiveness (Eq. [5]) and the constraints are dictated in view of the geo¬ 
metrical configuration (Fig. [2]): 


maximize ef[a,WN] 

{cx,Wm) 

subject to the constraints 

Ylj=o ~ ^j\ ~ ^ 

Ef=o«. = o 

lm[zi] =yi>{H- Hb) ( 8 ) 

~ ^ 

Im[z^] =y^=0, 

where the perimeter P, the height H and the height of the base Hb are assigned a priori. 
Note that the period is equal to one as it is used for non-dimensionalization, and the third 
constraint is an explicit equation to obtain the real constant R (Eq. [2]) . The fourth constraint 
defines the length of the fin and hence, the thickness of the plate. The equality Hf = H — Hb 
is achieved for sufficiently large perimeter P. 
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Similar to our previous work [32l |33l [Ml US HI [35] the integral Zi (Eq. [T|) is calculated 
using Gauss-Jacobi quadrature m, and we choose the WjS to be equispaced between 0 and 
1. The infinite product appearing in integral ([T|) can be truncated to a small value without 
affecting the accuracy due to the exponential decay of the hyperbolic sines m- The nu¬ 
merical optimization has been performed using the NLPQL optimization code developed by 
Schittkowski m- 

In Fig. H|we show numerical results of the fin effectiveness maximization problem (Eq. [8|) 
associated with a slab of height H = 0.5 and Hfy = 0.1, for different values of the perimeter 
(P) and Biot number (Bi). The results suggest that in all cases the optimum fin is infinitely 
thin, elongated in the vertical direction. It is very important to note that for a large Biot 
number the presence of the fin might not enhance the heat transfer rate from the base, i.e. 
the addition of a fin reduces the heat transfer rate. This can be justified as in the case of 
isothermal conditions the addition of a fin would result in a reduction of heat transfer (see 
Introduction, §1 first paragraph). 

The above results/conclusions suggest that: (i) the optimum fin is infinitely thin and 
elongated, and (ii) there exists a critical Biot number which characterizes whether a fin is 
effective or not. Hence, in the next section we determine the critical Biot number of a 
rectangular fin. 

3.1 Critical thickness characterizing the effectiveness of a fin 

The results obtained through the Shape Optimization procedure, suggest that for a fixed 
perimeter the optimum fin is an infinitely thin extended surface. It is tempting to infer that 
this is an artifact of the uniform heat transfer coefficient; had conjugate heat transfer been 
considered, the optimization procedure would have led to a different result. Furthermore, an 
interesting result that has emerged from the optimization analysis is that there exists a critical 
Biot number, associated with a given fin geometry, that determines whether the addition of 
the fin would enhance the heat transfer rate. In this section, we elucidate this point by 
considering the heat transfer rate associated with particular geometries, i.e. we consider the 
effectiveness of rectangular fin [2]. 

We consider the classical configuration of a periodic array of rectangular fins attached to 
an infinite rectangular base Fig [5k. Similar to the previous section, we non-dimensionalize 
lengths with the period of the array (L), to obtain the dimensionless Biot number Bi = hL/K. 
In Fig. [Sfi, we show a plot of the fin effectiveness Vs the Biot number for a rectangular fin of 
thickness Ht = O.I, attached to base of thickness Hi, = 0.1. The critical Biot number, i.e. the 
value of the Biot number where the effectiveness of the fin is equal to one, is BicriUcai ~ 16.4 
and it is independent of the length of the hn Hf. We have considered different sizes of 
rectangular hns and concluded that: 

1. A rectangular fin is effective if Bi Ht < 1.64. Expressed in dimensional vari¬ 
ables, fin thickness < 1.64 k/h. This is independent of both the thickness of the base 
and the length of the hn. 

2. The maximum effectiveness is realized at Bi = 0 and is equal to e/ [Bi = 0] = 

2Hj + l. 

The results have been verihed through hnite difference and hnite element simulations. 


(a) Bi = 1 



0 0.5 1 1.5 2 


X 


(c) Bi = 100 



0 0.5 1 1.5 2 


X 

Figure 4; Optimum surfaces/fins that maximize the heat transfer rate in a finite channel of 
height H = 0.5. The fins are assumed to be attached to a base whose minimum thickness is 
Hf = H - Hb = 0.1. In Fig. (a) we show results for Bi = 1 and two different perimeters: 
P = 1.24 (dashed curve, £f = 1.0); P = 1.68 (solid curve, £f = 1.4). In Fig. (b) we show 
results for Bi = 10 and two different perimeters: P = 1.15 (dashed curve, £f = 0.5); P = 1.63 
(solid curve, e/ = 1.0). In Fig. (c) we show results for Bi = 100 and two different perimeters: 
P = 1.27 (dashed curve, e/ = 0.5); P = 1.72 (solid curve, e/ = 0.99). 
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Figure 5; Fig (a): A periodic array of rectangular fins. Fig (b): Fin effectiveness (ej) Vs Biot 
number {Bi) for a rectangular fin of Hi, = Ht = 0.1. The three curves correspond to three 
different lengths {Hf). The lower curve corresponds to Hf = 0.01 while the middle and top 
curves to Hf = 0.1 and Hf = 1, respectively. 

4 Conclusions 

In this work we consider the heat transfer, shape optimization problem associated with a 
periodic array of extended surfaces subjected to convection with a uniform heat transfer 
coefficient. We address both the heat transfer problem and the Shape Optimization, inverse 
design problem, of finding the geometry that maximizes the heat transfer rate. 

The problem is formulated as a two dimensional, arbitrary channel of unit length which is 
bounded from above by a flat isothermal surface and from below by periodic extensions/fins 
with a uniform heat transfer coefficient. Given a set of geometrical constraints that charac¬ 
terize the geometry of the fin and the base, the objective is to find the optimum shape of the 
hn such that the heat transfer rate is maximized. Within the approximation of uniform heat 
transfer coefficient, the optimization procedure suggests that the optimum fin is inhnitely thin 
and long. 

Furthermore, the optimization procedure has revealed a very interesting result. There is 
a critical Biot number that characterizes the hn effectiveness. For values of the Biot number 
less than the critical a hn enhances the heat transfer rate, while for higher values it attenuates 
the heat transfer rate, hence it is not effective. This is investigated further by considering the 
simplest case of a rectangular hn of uniform cross section. Numerical results have elucidated 
the hn effectiveness and are summarized as follows: 

1. A rectangular fin is effective if Bi Ht < 1.64. Expressed in dimensional vari¬ 
ables, fin thickness < 1.64 k/h. This is independent of both the thickness of the base 
and the length of the hn. 

2. The maximum effectiveness is realized at Bi = 0 and is equal to £f [Bi = 0] = 

2Hf + l. 
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